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A phenomenological a-cluster model based on the charge symmetry of nuclear force allows one 
to estimate the last proton position radius (LPPR) in a symmetrical nucleus. The values of LPPR 
obtained for the symmetrical nuclei with 5 > Z > 45 are used in a long standing problem of 
determination of the Woods-Saxon single particle bound state potential parameters. With respect 
to the charge symmetry of nuclear force a requirement of equality of the nuclear potentials for the last 
neutron and the last proton in a symmetrical nucleus is added to the standard well-depth procedure 
in solving the Shrodinger equation for the nucleon bound states, which makes the Coulomb radius 
the crucial parameter to determine the others and the value of the last proton rms radius. The 
Coulomb radii have been obtained with using LPPR. Analysis of the last proton rms radii, the 
Coulomb radii, the nuclear potential radii obtained at the calculations in comparison with LPPRs, 
the Coulomb radii from the a-cluster model together with the experimental radii shows that for the 
nuclei with Z >15 it is inappropriate to represent a single particle bound state by the Woods-Saxon 
potential. For the nuclei with 5< Z <14 the error of the spectroscopic factor obtained with standard 
parameters in DWBA analysis of pure peripheral one nucleon transfer reactions is estimated. It is 
shown that for some nuclei using the standard parameters brings an error more than 20%. 

PACS numbers: 21.60.-n, 21.60.Gx, 24.50.+g 



I. INTRODUCTION 

It is well known that nuclear properties are described 
in the framework of three different representations of nu- 
cleon arrangements [l[ , the mean field representation, the 
liquid drop model and the a-cluster model. In anal- 
ysis of one nucleon transfer reactions A(a, b)B where 
A = B + n/p and b = a + n/p, one has to use the mean 
filed representation where neutrons and protons move in- 
dependently in the mean field. The nuclear mean field 
is described by the Woods-Saxon (W-S) potential due to 
the single particle model [2| with parameters of radius 
ro , diffuseness ciq and depth Vq depicting experimen- 
tal charge distribution. The standard Distorted Wave 
Born Approximation (DWBA) approach [1, Q , used in 
the analysis of differential cross sections of one nucleon 
transfer reactions, suggests that the nucleon wave func- 
tion with a factor S 1 / 2 where S is the spectroscopic factor 
(SF), is equal to the overlap function of the nuclei A and 
B. S is obtained by dividing the experimental cross sec- 
tion of the reaction at the main scattering peak by the 
calculated one. The experimental binding energy of the 
nucleon (neutron/proton) e„/ p in a nucleus A is used as 
the eigenvalue in the Shrodinger equation for the bound 
state. 

The determination of the single particle bound state 
potential parameters is a long standing problem (5rJl5j|. 
As the main criteria for selecting proper parameters the 
experimental nuclear radius and charge distribution are 
usually used @, UHl] with strong needs of additional 



model parameters like sizes of internal shells and the 
binding energies of the nucleons of the shells, which do 
not make the obtained values valuable. During a few 
decades in the DWBA analysis so called standard pa- 
rameters have been used, these are around the values 
ro = 1.25 fm, clq = 0.65 fm, which provide close values of 
the neutron and proton depths Vb w 50 MeV in so called 
well de pth procedure adjusting the depth to the values 
£ n/p 0tII|- Some articles, for example [13, EH, a re de- 
voted to finding so called global parameters of the W-S 
potential to describe single particle bound states of the 
nuclei united on some signs like closeness of the numbers 
of neutrons and protons to the magic ones. The parame- 
ters turned out to be close to the standard ones and they 
provide description of single particle spectra within an 
accuracy of 1-2 MeV [HI]. In this connection the task of 
finding the realistic values of the parameters which are 
supposed to be unique for every nucleus as well as an 
estimation of the error in obtaining SF by means DWBA 
with standard parameters is very actual. 

In [l2j the task of finding the proper values of the 
parameters was redefined for the task of finding proper 
value of asymptotic coefficient (AC) b n / p of the bound 
state wave function on the basis of the approach [l6| . 
where the asymptotic normalization coefficients of the 
bound state (ANC) C^/ p was introduced in the analy- 
sis. According to the approach the cross section of pure 
peripheral reaction at the main scattering peak of the 
angular distribution contains an asymptotic part of the 
nucleon's wave function which for the last neutron and 
the last proton are [H], [H| 

iPnWr-foo = b n i l kh l l:j (ikr) (1) 
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where h\j(ikr) is the spherical Hankel function of the 
first kind, W_„ y(2fer) is the Whittaker function, n is the 
Coulomb parameter for the bound state, k = yJ2ue n / p , 
a is the reduced mass for the bound system A = B + n/p. 
AC b n / p is determined by the geometrical parameters of 
the potential. Then ANC C 2 /p is related with SF and 



AC as follows [H 



r 2 

^n/p 
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This approach allowes one to obtain ANCs from anal- 
ysis of experimental cross sections of pure pe ripheral re- 
actions like A(d, t)B and B( 3 Re, d)A\XN$. The value 
ANC obtained from the reactions does not depend of the 
potential parameters used in calculation, because in that 
case S is in inverse proportion to the b^,. In [l2[ it 
was proposed to obtain AC and corresponding parame- 
ters from analysis of experimental cross sections of direct 
reactions with using a known value ANC. But later it was 
shown 17] that peripheral reactions are not good for the 
task. In case of reactions like (p, d) and (d, n) there is a 
strong dependence of the calculated cross sections on the 
optical potentials used for input and outgoing channels, 
which does not allow one to solve the task either. 

In (13| a method was developed where the parameters 
are obtained under the condition of the exact equivalence 
of the neutron and proton potentials for symmetrical and 
mirror nuclei (EPN condition) with respect to the charge 
symmetry of nuclear force. Together with the well- depth 
standard procedure the EPN condition makes the param- 
eters dependent on each other, so the Coulomb radius 
r DWBA k eC onres the main critical parameter to dctcr- 
min the value of root mean square (rms) radius calcu- 
lated for the last proton < r 2 > x / 2 and AC. Varying one 
of the parameters, for example diffuseness do, at r D ] WBA 
fixed brings some particular values of ro and Vq so that 
< r 2 > x l 2 ch ang es within 1% and b n / p changes within a 
few per cent Moreover, the ratio of squared neutron 
and proton ACs for mirror and symmetrical bound states 
b^/bp is quite a stable value (with variation within 2%) at 
a wide variation of the Coulomb radius [TEj . Some expla- 
nation why the ratio is to be stable is offered in [2(| ■ The 
value b 2 n jb 2 p can be used to predict the ratio C^/C 2 for 
mirror and symmetrical bound states, because the SFs in 
the framework of the shell model [2l[ are equal. In [lij to 
obtain some criterion to determine the Coulomb radius 
for the nuclei of lp shell some assumptions were made on 
the size of Is shell, which diminished the validity of the 
obtained values. 

In the present work like in [l1i - fl5| the EPN condition is 
used to obtain the potential parameters and AC. Besides 
£ n /p, which is calculated as the difference between the 
binding energies of the nuclei A and B, another macro- 
scopic quantity is used. That is the last proton position 
radius (LPPR) R p in the center of mass system (cms) of 
the symmetric nucleus A = B + n/p. The value of LPPR 
is estimated from analysis of the nuclear binding energies 
in the framework of an a-cluster model based on np-pair 



interactions with using charge symmetry of nuclear force 
acting between the nucleons belonging to the a- cluster 
matter f22T - |27j . One of the main aims of developing the 
model was finding proofs of validity of the EPN condi- 
tion for single particle bound states in symmetrical and 
mirror nuclei. 

The Coulomb radius of the last proton interaction with 
the residual nucleus R% p VBA = rg WBA (A - l) 1 / 3 is de- 
fined from the formula used for the Coulomb potential in 
the standard DWBA (see the manual to DWUCK pro- 
gram (U) 



{Z-l)e 2 

OUDWBA 
2K Cp 
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at r < Rg p VBA and 



(5) 



at r > R D ] p VBA where the Coulomb energy E p of the last 



proton 



AE„ 



(6) 



and AE np = e n — e p . Eq. (6) is the consequence of 
the EPN condition. Taking into account the observed 
fact that the Shrodinger equation at the EPN condition 
provides the value of rms radius < r 2 > x / 2 less than 
r dwba Eq ( 4 ) ig used tQ find the value R dwba ( Qr 

DWBA\ 



P 

the parameter Tq WBA ) at r = R p 



The article consists of 4 sections. The second section 
explaines how the values LPPR R p are obtained in the 
framework of the a-cluster model for the symmetrical nu- 
clei with 5< Z <45. In the third section the parameters 
obtained by means of the EPN method are given. In the 
section 4 there are discussions of the results. 



RLPP OBTAIND IN ALPHA-CLUSTER 
MODEL 



II. 



The a-cluster model 22 i27j is based on the charge 
symmetry of nuclear force, so the Coulomb energy of the 
last proton in a symmetrical nucleus is to be equal to 
AE np (6). According to representation of an a-cluster 
liquid drop the nuclear matter is incompressible, so the 
Coulomb energy of a nucleus equals the energy of the 
np-pairs consisting the a-clusters. Then for the nucleus 
with even Z , with the number of a-clusters N a = Z/2, 
and for the nucleus with odd Z\ = Z + 1, the number of 
a-clusters N a + 0.5, the empirical value of the Coulomb 
energy is supposed to be estimated as follows [25[ 



E l 
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where AE np i and AE np are the differences of one np- 
pair neutron and the proton binding energies, odd np- 
pairs have index 1. From the analysis of the experi- 
mental binding energies and the Coulomb energies (7) 
of the lightest nuclei with Z, Z\ < 8 some important val- 
ues were obtained [13, HH ■ These are the binding energy 
and the Coulomb energy of one a-cluster (absolute val- 
ues are given) e a = £4 Hc (£4 Hc =28.296 MeV (28|) and 
£^=0.764 MeV, the binding energy and the Coulomb en- 
ergy of a cluster-cluster interaction e aa =2.425 MeV and 
e cL~ 1.925 MeV and the Coulomb energy of the inter- 
action of the single np-pair in odd Z\ nucleus with the 
a-cluster of its close vicinity s^ pa = 1.001 MeV. Also a 
simple formula to describe the binding energies of sym- 
metrical nuclei with Z, Z\ <29 was found 22] E = 
N a e a + 3(iV Q — 2)e aa for the nuclei with even Z and 
Ei = E + 14 MeV for the nuclei with odd Zi, which 
means that one added a-cluster brings three new links 
with the closest ones. The formula means that the long 
range Coulomb energy must be compensated by the sur- 
face tension energy, which allows one to find a formula 
to calculate the latter. Thus a successful formula to cal- 
culate binding energy for stable and beta-stable nuclei 
as well as for the nuclei around the stable valley has 
been found (2a |. In the formula the binding energy of 
a-clusters is calculated separately from the energy of ex- 
cess neutron pairs (nn-pairs). The accuracy of the cal- 
culation is comparable with the Weizsacker formula [2i| , 
but unlike this well known formula the parameters used 
in the model are not fitting ones, they had been found 
from analysis of the binding energies of reduced amount 
of symmetrical nuclei. Another important finding prov- 
ing the a-cluster representation is that the nuclear radius 
for a stable and a beta-stable nucleus is defined by the 
number of the a-clusters rather than by the total number 
of the nucleons. The simplest formula for nuclear radius 
growing with Z, Zi is as follows 0, HH| 



R = R a N 1 J 3 ;Ri = R a (N a + 0.5) 



1/3 



(8) 



where R a = i?4 Ho (i?4 Hc =1.71 fm [30] ) for the nuclei with 
Z=2 and 5< Z, Z x < 10. For the nuclei with Z, Z x > 24 
R a =1.595 fm. To widen the number of nuclei to be de- 
scribed to the ones with N > Z, Z\ the model was devel- 
oped to the representation of nucleus as a core (liquid a- 
cluster drop with dissolved nn-pairs in it) and a molecule 
of a few a-clusters on its surface [2(| [27| . The notion of 
a nuclear rnolecule on the surface of a core was devel- 
oped in [3l[. The number of the molecule's a-clusters 
is obtained from analysis of the nuclear binding energy. 
So for the nuclei with 10 < Z,Z\ < 24 in (8) the ra- 
dius R, Ri is defined by the sum of the volumes of the 
surface molecule, presumably the nucleus 20 Ne and 23 Na 
(in case of Z\ nuclei one excess neutron is glued to the 
single np-pair) , and of the growing core consisting of the 
a-clusters of the radius R a =1.595 fm [2(| . The surface 
tension seems to be responsible for existence of core. It 
was shown that for the most stable nuclei the specific 
density of the core binding energy p is an approximately 



constant value p — 2.55MeV/fm 3 at the number of the 
surface molecule a-clusters equal to three (three and a 
half with one excess neutron in case of odd Z\ nuclei). 
This provides an explanation of the particular number of 
excess neutrons in stable nuclei [27| . In the nuclei with 
smaller amount of excess neutrons the core is smaller (the 
molecule is bigger), because the number of neutrons pro- 
vides the p (a little bit less than the saturated value) for 
a smaller number of a-clusters. This theory allows one 
to calculate the radii for the nuclei with N > Z from 
analysis of the binding energy (26|. It explains the phe- 
nomena of a slight grow of the radii of the isotopes of 
one Z, Zi with decreasing A, A\. The obvious success of 
the model in describing binding energies and radii of the 
nuclei proves the validity of the EPN condition and its 
consequence (6). 

The simplest way to obtain the value of LPPR in the 
cms of the nuclei A = 2Z and A\ = 1Z\ is given by the 
formula 
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AE np i 



(Z-l)e 2 



A-l 



(Z, - l)e 2 



R 



Ai-l 
pi 
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where R p 1 is the radius of the last proton position in 
the cms of the residual nucleus with the mass A — 1. In 
case of Rp' 1 < R^" 1 where R^T 1 = r c (A - l) 1 / 3 the 
following equation is used (see (4)) 



AE, 
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where rc is obtained from the formula for the Coulomb 
energy Ec, Eci (7) of the charge sphere with the radius 
Rc = rcA 1 / 3 and R C i = r ci ^ /3 



E c = 3/5 
Ef = 3/5 



Z 2 e 2 

Zle 2 
Rci 



(11) 



Values R p and R p i are related with R£ 1 , R^ 1 as fol- 
lows 



Rp Rp 



a-i(A-I) 



A 



JXpi — n pl ^ 



(12) 



The values R P:P i obtained by means (9) and (10) are 
given m Table U For the nuclei with Z,Z X <8 the val- 
ues R PtP i obtained by (10) are considerably less than 
R exp , which is out of reason, and they are not presented 
in the table. (For example R p =2.209 fm for 16 O and 
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_Rie =2.718 fm [30]). For light nuclei the representation 
of nucleus as a charge sphere (11) is not good, therefore 
the values rc,d do not give reasonable results to be used 
in (10). 

The a-cluster model also allows one to estimate the 
value of LPPR upon a suggestion that the Coulomb en- 
ergy of the last proton (or np-pair) in even Z nucleus 
comes from the Coulomb interaction between two protons 
(two np-pairs) in the last a-cluster P ms the Coulomb 



TABLE I: The RLPP R PtP i calculated in the framework of 



a-cluster model representation. E 
(7), rc is the Coulomb radius (11) 



is the Coulomb energ 

AEnr> — — En S"p [28 



The values R p 9 pl and R pp i are obtained by means Eq.s (9) 

and (10) correspondingly with using (12), the values arc 

calculated by (13) together with (15), R^pl are calculated by 
(14) with (16), Rp"p{ is calculated by (17). 
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.684 


2.763 


2, 


.389 


2 


.491 






2 


.491 






7 


11.070 


1 


.587 


3.003 


2. 


.672 


2 
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14.607 
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.502 


3.537 


2. 


.672 
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.764 
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18.149 
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.471 


3.542 


3. 


.072 


3 


.039 
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.932 
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.919 


3. 


266 


10 


22.170 


1 


.436 


4.021 


3. 


062 


3 


.102 


2 


.872 
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.172 


3. 


263 


11 


26.500 
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.408 


4.330 


3. 


.174 
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.117 
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.021 
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.050 
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279 


12 


31.338 
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.376 


4.838 
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.138 
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.133 
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.927 
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.231 
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13 


36.397 
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.354 


5.059 


3. 
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.234 
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.143 
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.181 


3. 


.378 


14 


41.992 
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.328 


5.595 


3. 
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.219 
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.021 


3 


.310 
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.327 
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47.717 
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.311 
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.369 
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.286 
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53.896 
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.293 


6.179 


3. 


.386 
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.397 
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.236 
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.471 


3. 


.483 


17 


60.261 
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.279 


6.365 


3. 


.513 
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.490 


3 


.410 


3 


.441 


3. 


.620 


18 


67.010 


1 


.265 


6.749 


3. 


.526 


3 


.547 


3 


.410 


3 


.612 


3. 


.621 


19 


73.940 


1 


.255 


6.930 


3. 


.642 


3 


.635 


3 


.559 


3 


.581 


3. 


.764 


20 


81.256 


1 


.244 


7.316 


3. 


.646 


3 


.680 


3 


.548 


3 


.731 


3. 


.759 


21 


88.533 


1 


.238 


7.277 


3. 


.863 


3 


.893 


3 


.825 


3 


.820 


4. 


.035 


22 


96.182 


1 


.232 


7.649 


3. 


.864 


3 


.934 


3 


.812 


3 


.958 


4. 


031 


23 


04.090 


1 


.225 


7.908 


3. 


.919 


3 


.944 


3 


.870 


3 


.881 


4. 


081 


24 


12.316 


1 


.219 


8.226 


3. 


.942 
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.005 
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.888 
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.034 
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25 
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.213 
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.953 
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.207 
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.012 
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.071 
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.954 
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.101 
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38.698 
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.201 
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.046 


4. 


.069 
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.987 


4 


.018 


4. 


.202 


28 


48.173 
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.195 


9.475 
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.030 
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.085 
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.955 
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.113 


4. 


.187 


29 


57.727 
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.190 


9.554 


4. 


.147 


4. 


.180 
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.099 


4 


.127 


4. 


.314 


30 


67.606 


1 


.185 


9.879 


4. 


.157 


4. 


.220 


4. 


.101 


4 


.243 


4. 


316 


31 


77.643 


1 


.181 


10.037 


4. 


.235 


4. 


.273 


4. 


.192 


4 


.220 


4. 


407 


32 


88.033 


1 


.176 


10.390 


4. 


.229 


4. 


.297 


4. 


.175 


4 


.314 


4. 


401 


33 


98.243 


1 


.174 


10.210 


4. 


.445 


4. 


.511 


4 


.430 


4 


.443 


4. 


.660 


34 


09.196 


1 


.170 


10.953 


4. 


.275 


4. 


.371 


4. 


.213 


4 


.357 


4. 


544 


35 


19.807 


1 


.168 


10.611 


4. 


.548 


4. 


.621 


4. 


.538 


4 


.552 


4. 


.774 


36 


30.729 


1 


.167 


10.922 


4. 


.550 


4. 


.653 


4. 


.535 


4 


.647 


4. 


.777 


37 


41.984 


1 


.164 


11.255 


4. 


.544 


4. 


.606 


4. 


.521 


4 


.546 


4. 
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38 


53.367 


1 


.163 


11.383 


4. 


.619 


4. 


.706 


4. 


.603 


4 


.715 


4. 


.799 


39 


65.065 


1 


.160 


11.698 


4. 


.618 


4. 


.683 


4. 


.596 


4 


.624 


4. 


.828 


40 


76.819 


1 


.159 


11.754 


4. 


.718 


4. 


.805 


4. 


.706 


4 


.816 


4. 


.894 


41 
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1 


.158 


11.830 


4. 


.810 


4. 


.895 


4. 


.803 


4 


.827 


5. 


.056 


42 


300.617 


1 


.158 


11.968 


4. 


.874 


4. 


.983 


4. 


.870 


4 


.978 


5. 


102 


43 


312.857 


1 


.157 


12.240 


4. 


.884 


4. 


.972 


4. 


.877 


4 


.904 


5. 


.134 


44 


325.385 


1 


.156 


12.528 


4. 


.886 


5. 


.000 


4. 


.877 


4 


.985 


5. 


139 


45 


338.023 


1 


.155 


12.638 


4. 


.958 


5 


.048 


4. 


.952 


4 


.981 


5. 


.211 



long range interaction of the last np-pair with the resid- 
ual nucleus of the mass number A — 4, consisting of N a -1 
a-clusters. For odd Z\ nucleus the energy of the last 
proton comes from the energy of the last np-pair interac- 
tion with the closest a-cluster e^ pa plus the energy of its 
Coulomb long range interaction with the residual nucleus 
of the mass number A-6 [251 ] 
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(13) 



pi 



are the distances between the 



last np-pair and the cms of the residual nucleus with 
mass number A — 4 and A\ — 6. 

Another formula comes from the a-cluster model rep- 
resentation that the last a-cluster has three links with 
the nearest a-clusters. Then the last a-cluster Coulomb 
energy AE a is equal to the sum of its own Coulomb en- 
ergy e^, the energy of the three links with the nearest 
clusters 3e^ a and the energy of the long range interaction 
with the rest a-clusters of the nucleus 2(Z — 8)e 2 / R p ~ 16 , 
where R p ~ 16 is the distance between the last proton and 
the cms of the remote a-clusters with total mass A — 16. 
In case of odd Z\ nucleus the Coulomb energy of the last 
np-pair is the sum of the Coulomb energy of its inter- 
action with the three nearest a-clusters 3e^ pQ and the 
Coulomb energy of the long range interaction with the 
rest a-clusters of the nucleus (Zi — 7)e 2 / R p l~ 1A , where 
i?p 1-14 is the distance between the last proton and the 
cms of the remote a-clusters with the total mass number 

A x - 14 m 



AE a =s% + 3e, 



.c 



AE. 



npl 



3e 



2(Z~8)e 2 

rip 
{Zi - 7)e 2 



npa 



R 



pi 



(14) 



where AE a = AE np x + AE np . The simplest formula for 
the Coulomb energy decreasing with distance is used, as 



the values R 



p A - 4 , Rpx 1 " 6 : Rp~ W and i?^ 1 " 14 are surely 
bigger than the nuclear radius. 

The distance between the cms of the last a-cluster and 
the cms of the residual nucleus with mass number ^4-4 is 
estimated as i?^" 4 - i?.4 Ho . Then (R A ~ 4 - R4 Hc ){A - 
4) I A is the radius of the last a-cluster position in the 
cms of the nucleus A. The same logic is used for the 
odd Z-i nucleus. (R Al ~ 6 - i?4 HG +np)(^i - fy/Ai is the 
radius of the position of the last cluster consisting of the 
1.5 a-clusters in the cms of nucleus A\ and i?4fj c +np = 
i?4 Ho 1.5 1/3 (8). So LPPRs in cms of nucleus A R PiP i are 
calculated as follows 



R p = (R p 



A-4 



A - 4 

-R 4 Hc) -, 1" R^Hc 



A 



R 



pi 



( R pl 



Ai-6 
pi 



R- 



He+np ) 



A x -6 
~Ax~ 



R 



Hc+np- 



(15) 
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To calculate Rp, p i in case of (14) one can use LPPR in the 
nuclei 16 Rf and 14 N obtained by (9), see TableQ] 
Then using the same logic as in (15) we have 



Rp — (R, 



A-16 
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R p ) 



p _ /dAi-14 
ri pl — (K pl 



R 



pi) 



A- 16 
A 

Ai - 14 
~A^~ 



Hi 



R 



pl4 

n pl . 



(16) 



latter is that R^li are obtained without any parameters 
only upon the suggestion (6). For further calculations 

(9) 

the values LPPR R PtPl are used and they are indicated 



as R 



p. pi ■ 



III. BOUND STATE POTENTIAL 
PARAMETERS IN SYMMETRICAL NUCLEUS 



Average value Rp V p[ of the radii estimated with using 
the a-cluster model parameters rc, £„, e^ a , e% pa is also 
calculated as follows 



reaver 
^p,pl 



/d(10) 



d(13) 
^p t pl 



The values LPPR in cmc of nuclei A, A\ estimated by 
different ways are given with upper indexes correspond- 
ing to the Eqs, see Table I. One can see from the table 
that the values R p ^ v \ obtained without model parame- 
ters but with the assumption (6), which is the straight 
consequence of the EPN condition, are in agreement 
with Rp V p[ with the average deviation in 0.065 fm. For 
some of the light and the heaviest nuclei the difference 



\R 



(9) 



R° 



is more than 0.1 fm. 



Let us check if the values R^i and Rp V p[ are consis- 
tent with the experimental radii. Not for all symmetrical 
nuclei there are experimental data. In those cases the ex- 
perimental radius of the nearest isotope is used [32[ . In 
the case of no data, this is the case of Zi=43, Eq. (8) is 
used. The squared nuclear radius ((R A ' Z ) 2 , (R Al ' Zl ) 2 ) is 
calculated as the sum of the squared radius of the residual 
nucleus ((i?- 4 " 1 ^ 1 ) 2 , (R^- 1 ' 21 ' 1 ) 2 ) and the square ra- 
dius of the last proton position in the cms of the residual 
nucleus ((i?^" 1 ) 2 , (-R^ 1 _1 ) 2 ) weighed 



{R A, Z) 2 



Z 



1 



(Zi 1) 
Zi 



1 



(i^-^-y + ^(< i - l ) 2 <i8) 



av p e {\ does not 



With growing Z, Z 1 the difference \R p ^ pl — Rp 
affect much the calculated nuclear radii (18). The aver- 
age deviation \R exp — R t h\ where R t h = R A ' Z , R Al ~ 



(9) 

p,pi 



Hh 

and 0.056 fm for R%™[. For 



is 0.058 fm for R, 

the nuclei with Z,Z\ = 5,6,7,8,9,10,11 the difference 



\R 



(9) 
■p,pl 



R 



p,pl 



0.272, 0.102, 0.189, 0.092, 0.033, 0.040, 



0.057 fm leads to the corresponding differences in the val- 



ues \R e 



R 



th 



calculated with using R pp x and R p ^ p l 



0.069, 0.002, 0.032, 0.001, 0.005, 0.005, 0.006 fm, which 

l pl | give close values of 



shows that using R p 9 pl and i?™ er 



Rth- It should be noticed here that for the nuclei with 
Z, Zi =5, 6, 8, 10, 11 R p V p{ gives better agreement with 
the experimental data. For the nuclei with Z, Z\ < 8 
R { (see Table I) . It is obvious that 



the values R£™{ 



(13) 

in the case of the light nuclei the values R p p i are more 
preferable than R p li , but the obvious advantage of the 



Having obtained LPPR R p . p i one can numerically cal- 
culate the radius R Cp = rg WBA (A - l) 1 / 3 and R Cp i = 
r DWBA^ Ai _ ^1/3 from ( 4 ) at E(jp = AEnp and r = 

(17) Rp.pi- The obtained values r^^ BA are used to solve the 
Shrodinger equation for one nucleon bound state with the 
W-S potential with the spin-orbit part of Thomas form 
with A=25 as it is used in DWUCK program Q 



V(r) = V (f(r) + 



A ldf(r) 



45.2 r dr 



La) 



where 



f(r) = [1 + exp(- 



'Y) 



(19) 



(20) 



The EPN-condition used together with the well-depth 
procedure makes the parameter r^^, BA the critical one 
in determination of the last proton rms radius < r p pl 



> 



and AC 6„ jP [13J, |14||. Variation one of the parameters at 
the r^Qi fixed, for example the diffuseness ao = 0.4-5- 
0.7 fm, changes the other parameters ro and Vq that way 
that < r 2 ,„ > x / 2 stays almost the same (variety within 



%/p 



1%) and AC b n / p changes within a few per cent [14|. So 
one can say that in such calculations the Coulomb radius 
defines last proton rms radius. Quantum numbers n, l,j 
are selected according to the Pauli conservation principle, 
parity conservation rule and the sum rule for momenta. 
In Table II for the cases when both binding energies and 
spins of nuclei A and B are known [28| , the values 



„DWBA 
C.Cl 



and corresponding values ro and Vq are given together 
with the rms radii for neutron/proton < r 2 y p >*/ 2 and 
AC b n / p . 

Here ao=0.65 fm is used. The value ao=0.55 fm in 
the case of nucleus 24 Mg gives r =1.2592 fm Vb=60.0168 
MeV < r 2 > 1 / 2 =3.160fm, < r 2 > 1 / 2 =3.210 fm, 
fr n =8.561/m -1 / 2 and 6 p =12.54 /m -1 / 2 . Comparison of 
the values with the data in Table II shows that ao=0.55 
fm almost does not change rms radius and AC changes 
for 6%. Searching parameters under the EPN condition 
is done by a special program with using iteration. Re- 
sult may depend on the starting values, which shows the 
accuracy of the calculation. This is few 0.001 fm for 
r and few 0.01 MeV for Vq, which consequences the 
rms radius accuracy of 0.001 fm and the AC accuracy 
of few 0.01 fm -1 / 2 . For example, in case of the nucleus 
Mg the result does not depend on the starting values. 
For the nucleus 44 Ti another solution gives ro=1.2719 
fm and Vo=-60.8539 MeV, < r 2 > 4 / 2 =3.962 fm and 
bp = SlAtfm- 1 / 2 . 



6 



TABLE II: Parameters of the W-S potential at the EPN-condition. In the first column the bound state A — B + n/p, in the 
2nd and the 3d columns there are the quantum numbers n,l,j and experimental binding energies e n ,j>, the 4th column gives 
Coulomb radii r^^Q^ A (4), the 5th and the 6th contains the parameters ro and Vb at do =0.65 fm and rms radii < r\ „ > 1,/2 , the 
7th column contains -R P , P i (see Tabled, the 8th column presents AC b n , p , the 9th contains AC b„ iP calculated at the standard 
parameters, the 10th column gives the error brought by the standard parameters in DWBA analysis to obtain SF (21). 
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The procedure of seeking parameters under the EPN- 
condition brings the determination of rms radius < 
r p,pi > an d radius r by the Coulomb radius r BV ^, BA 
in a way that decreasing the latter causes increasing the 
values < r ppl > 1 / 2 and r . 

In the table also the values AC b^ p are given which 
have been calculated at the standard parameters ro=1.25 
fm, a =0.65 fm with rg^f A =1.25 fm at well-depth pro- 
cedure without addition of the EPN condition. One 
wants to know what error in obtaining SF in analysis 
of one nucleon transfer reactions by means of DWBA 
is brought with using standard parameters. In case of 
the pure peripheral reactions the cross sections at main 
scattering peak contains an asymptotic part of the wave 
function (1) and (2). Then the deviation of the SF S ex- 
tracted by the DWBA analysis with standard parameters 
from the value obtained with using the EPN-parameters 
is calculated as follows 



K 



err 



/P 



iK /P ) 2 \ 



* 100% 



(21) 



i/p 



The values R™p{ (see Table I) for the light nuclei have 
also been tried in obtaining the parameters. In this case 
the proton rms radii are also in agreement with the values 
within 0.1 fm. For the nuclei with Z, Z x = 5, 6, 7, 8, 
9, 10, 11, 12 the parameter rg^f^=1.955 fm, 1.360 fm, 
1.553 fm, 1.312 fm, 1.510 fm, 'l.389 fm, 1.434 fm, 1.335 
fm, which brings the r , V = 0.88 fm, -79 MeV; 1.24 fm, 
-64 MeV; 0.93 fm, -87 MeV; 1.38 fm, -50 Mev, 1.12 fm, - 
69 MeV; 1.28 fm, -51 Mev; 1.14 fm, -64 MeV; 1.23 fm,-64 
MeV. This means that the parameters given in Table II 
should be taken with the probability corresponding to the 
deviation between the results of using and R™ p [ , 

see Table III. It is seen from the table that deviation in 



TABLE III: Deviations in results of searching parameters with 
using Rp, p i and Rp^p{ ■ The indications are the same as in 
Table II.' 
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r c l ci A - Varying r^f A = l.l-1.5 fm brings a change of 
the ratio within 2 %, which can be used to predict the ra- 
tio of experimentally obtained values ANC's C^/Cp due 
to the equality of the spectroscopic factors of symmetri- 
cal bound sates. 



IV. DISCUSSION OF RESULTS 

In FIG. [1] the graphs of the values of different kinds 
of radii in dependence on Z are given (in this section for 
further convenience even Z and odd Z\ are indicated as 
Z, as well as the other values like A, R p , Rc and so on) 
[t] 

These are the Coulomb radius Rc = rcA 1 ^ (11) of 
the sphere having the Coulomb energy E c (7), the values 
of LPPR the values R2 WBA = r ° WBA A 1 ^ where 



„DWBA 



is smaller for bigger Z, Z\ , so is the deviation in nuclear potential determined by the parameter 



DWBA 
1 C,C1 

the error. In case of the nuclei 10 B, 12 C, 8 and 10 Ne 
the smallest possible error in obtaining SF with standard 
parameters in case of pure peripheral reactions is small, 
which is 5, 3, 7, 7. But if the proper parameters diffcr 
in full length of the possible deviation the error becomes 
significant, except for 10 Ne. 

In the calculations the value b^/b p obtained under the 
EPN condition is almost stable value to a wide variety of 



up, uiic miura ±c c — i c 

• c is obtained by (4) with using R p , rms radius 
for the last proton < r 2 p > x / 2 calculated with the wave 
function obtained with using the EPN parameters for 



„DWBA 



Also the values of the radius of the nuclear potential Rq = 
roA 1/>3 (for ro see Table II) and for the standard values 
ro = rQ WBA =1.25 fm in comparison with experimental 
nuclear radii R exp are given. The nuclear radii R( A ' Z ' 
are calculated by (18) where instead of R A ~ X the proton 
rms radii transferred into the cms of the nucleus A-l by 
(12) are used. 

Comparison of the values Rc and R C >WBA shows that 
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Z -| i 1 1 1 , 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 . 1 , 

4 8 12 16 20 24 28 32 36 40 44 48 Z 

FIG. 1: Nuclear radii. Symbols indicate the following values: 
the experimental nuclear radii R eX p (pentagons) [281 ] (in case 
of absent data for the symmetrical nuclei the radius of the 
nearest isotope is taken, for the nucleus with Z=43, which 
is not presented in [2f|, the radius is calculated by (8)), the 
Coulomb radius Rc (11) (triangles), the radius Rq WBA — 
r DWBA^i/3 ^ (squares) and the radius of nuclear potential 
Ro ~ roA 1 ' 3 satisfying the EPN-condition (crosses). Lines 
indicate: RLPP R v (solid thick line), rms radii < r 2 > 1//2 
(thin line), calculated nuclear radius R A ' Z (18) (solid thin line 
gong along the symbols of R e x P )(in Eq. (18) instead R A_1 the 
values rms radii < r 2 > 1//2 transferred into the c.m.s. of the 
nucleus A-l by (12) are used). Also the standard potential 
radius Ro = 1.25 A 1 * 13 (dashed line) and the corresponding 
rms radii < r 2 > x / 2 (line of dots) are given. 



they are in agreement for the nuclei with Z=5 and 
9< Z <32 with the average deviation in 0.12 fm. The 
reason of the considerable difference between the values 
for the light nuclei Z < 9 is that the approach of a charge 
sphere is not good for the case. That is one of the reasons 
why the value rc is not used here to obtain the parame- 
ters. The other reason is that for the nuclei with Z > 32 
the values rc < rc WBA , so the corresponding potential 
has either too big radius or in some cases there are no 
parameters satisfying the EPN-condition. 

In FIG. [5] the Coulomb energy E c (7) of the charge 
sphere of radius Rc is given in comparison with the 
energies of two charge spheres with the radii equal to 
Rg WBA (R p ) and R° WBA {R% ver ) obtained by (4) with 
using values R p and Rp Ver . For the nuclei with 5< Z <26 
all three lines deviate within few MeV. For heavier nu- 
clei they fork. The Coulomb energy of the sphere of 
the radius R B ' WBA (Rp lver ) is the highest graph and it 
splits from the others at Z >26, and these two deviate at 
Z >32. E c is the middle line between the Coulomb en- 
ergies of the two radii spheres. For the heavier nuclei the 
deviation grows within a corridor E ±15 . For compari- 
son the Coulomb energy according to the Weizsacker for- 
mula [29| is given too, which is in agreement with the first 
three graphs within average deviation in few MeV for the 
nuclei 5< Z <32 and for heavier nuclei it is in agreement 
with the energy of the sphere of radius R c rwBA {R p ). 



350- 




4 8 12 16 20 24 28 32 36 40 44 48 Z, 



FIG. 2: The Coulomb energies of two charge spheres having 
the radii Rq WBA calculated by (4) with R p (lower line) and 

R DWBA = r DWBA A l/3 calculated with R aver ( upper ^ ^ 

comparison with the Coulomb energy Ec (7) (middle line). 
The Coulomb energy due to the Weizsacker formula [2^] is 
indicated with dotted line. 



The mean field model considers a nucleus as consists 
of neutrons and protons moving independently in a mean 
filed. The question about consistency of the W-S poten- 
tial used as a single particle potential with experimental 
data on nuclear charge and matter distributions was dis- 
cussed in a number of articles, for example 0,0-0). There 
was shown that theoretical charge distribution calculated 
with using the sum of squared single particle wave func- 
tions produced by the W-S potential with a radius close 
to the standard one is in agreement with experimental 
charge density distribution for the nuclei with large devi- 
ation of A. But the sizes of internal shells are not known, 
so the fitting experimental radius can be done by varying 
sizes of internal shells. Besides, as it is already shown in 
section II experimental radius alone can't be a sensitive 
criterion for checking validity of the last proton rms radii 
calculated at the parameters. 

The values LPPR obtained from independent analysis 
of nuclear binding energies allow one to check the rms 
radius < r 2 y 1 / 2 for the last proton produced by the 
W-S potential. From FIG. [T] (see also Table II) one can 
see that for the nuclei with Z < 14 the EPN potential 
provides a good description R p , the difference is ~ 0.1 
fm. For the other nuclei with Z > 14 the difference 
< > x / 2 —R p is larger, more than 0.2 fm. It leads one 
to a conclusion that the W-S potential as it is used in the 
single particle potential model is not good for the heavier 
nuclei. The values LPPR have also been tried here to 
obtain the Coulomb radius r ^ WBA from requirement of 
equivalence of R p —< ri > 1 / 2 (first it was proposed in 
[la ) for some of the nuclei of Id shell but was not fulfilled 
correctly, because the values of LPPR were taken in the 
cms of the residual nucleus). In that case for heavier 
nuclei the value rc WBA becomes unreasonably big or 
there is no solution of the Shrodinger equation at all, 
which leads to the same conclusion. 
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In spite of the fact that the proton rms radii for the 
nuclei with Z > 14 considerably deviate from the R p 
the nuclear radii RS A > Z ' (18) with using < r p > 1 / 2 (re- 
calculated (12) for the cms of the residual nucleus with 
mass number A — 1) are consistent with R exp , see FIGfT] 
That is because the relative weight of the last proton de- 
creases with growing Z. This is another evidence that 
the experimental nuclear radius alone can not be used as 
a sensitive test for validity of the last proton potential 
parameters without an additional criterion like R p . 

Another proof of unavailability of the W-S potential is 
the size of the nuclear potential well. The potential ra- 
dius Rq grows with Z by a periodic function, see FIG. [T] 
The values Rq become too big in comparison with the ex- 
perimental radii. For heavier nuclei the difference grows 
up to 2 fm, which is not right taking into account the 
short range nuclear force. For some nuclei the radius de- 
scends to the standard value 1.25A 1 / 3 , which is also too 
big in comparison with the experimental radius (the ex- 
planation is in next paragraph). One also should notice 
here that for the nuclei with Z > 15 radius Rq increases 
by a leap. Searching parameters under the EPN condi- 
tion for the nucleus 30 P gives V = 52.78 at i? =4.07 fm 
at tq — 1.31 fm, see Table II. Thus, one can suppose that 
the number of nuclei for which the W-S potential is good 
is restricted by those with Z < 14. For the nuclei with 
Z > 15 the radius Rq obtained at the EPN condition 
growing with a leap reveals that the W-S potential is not 
proper in the framework of single particle representation. 

Let us consider the consistency of the W-S potential 
with experimental charge radii at the standard param- 
eters used without the EPN condition. The difference 
between LPPR R p and < r 2 > 1 / 2 at the standard pa- 
rameters for the nuclei with Z < 14 is bigger than that 
for the parameters of the EPN condition, see FIG. [U 
For the odd Z nuclei it exceeds 0.2 fm. For the nuclei 
with Z > 14 the rms radii are in the better agreement 
with LPPR. But then the potential radius 1.25 A 1 / 3 is 
not consistent with the experimental radii. Experimen- 
tal radius 1.005 A 1 / 3 < R exp < 1.077A 1 / 3 for the nuclei 
with 10< Z <24 and R exp = 1.005 A 1 / 3 for the other nu- 
clei with Z > 24, see Eq. (8) rewritten for A = 4N a . It 
means that the nuclear potential radius 1.25A 1 / 3 grows 
with Z with a bigger rate than the radius of the nucleus, 
which is in a clear discrepancy with the short range nu- 
clear force. One can see in FIG. [1] that the Coulomb 
radius Rc obtained in the framework of the a-cluster 
model relates with R exp as Rc = R exp + d where d is an 
approximately constant value. Potential radius Rq can 
be equal to Rc or, if it is different, at least it should 
grow with Z in consistency with Rc and with R exp . So 
one can suppose that the relation between Rq and R exp 
should be like Rq = R exp + d, which is consistent with 
the nature of short range nuclear force. 

Besides, at standard parameters the neutron and pro- 
ton potential wells are allowed to be different. At tq = 
1.25 fm in case of the nucleus 30 P neutron's potential 
depth is less than the proton's one, Vo n = 57.15 MeV 



and Vop = 57.50 . The difference does not seem signifi- 
cant in point of view of mean field theory where neutrons 
and protons are distributed independently and their cen- 
ters of mass do not coincide. But it is not right in the 
representation of the a-cluster model, where a nucleus 
consists of np-pairs joined in a-clusters and the posi- 
tion of the last neutron and the last proton belonging 
to one pair is determined by one potential. Calculations 
show that at the EPN condition the proton stays farther 
from the center of mass of a nucleus A than the neu- 
tron < r 2 > 1/2 >< r 2 > 1/2 (see Table II), which also 
appears natural. If one takes into account the difference 
in neutron and proton masses, the neutron potential is 
supposed to be deeper than proton's one, but not the 
opposite. It should be pointed out that the standard 
potential radius 1.25A 1 / 3 is less than the potential ra- 
dius Rq obtained at the EPN condition only because at 
the standard parameters the neutron and proton depths 
are allowed to be different at unreasonable prevalence of 
proton's nuclear potential above neutron's one. Thus, 
one can see here that a good description of experimental 
charge radii in the framework of single particle poten- 
tial model is provided by an unreasonably big radius of 
the potential at both standard parameters and at the pa- 
rameters of the EPN condition. It is not so big as one 
at the EPN condition, because the neutron and proton 
potential wells arc allowed to be different in the standard 
well-depth procedure. 

There is another remark about not consistency of the 
W-S potential with the nuclear density distribution. In 
the self-consistent calculations it is shown [33| that the 
single particle potential should have no symmetry in 
the surface thickness to be consistent with the nuclear 
density. The internal part is to be considerable larger 
to.5 ~ ^0.9 > to.i ~ to.5 where to. 5 means the half-potential 
radius. 

If for the light nuclei with Z < 14 the W-S potential is 
good, one has an opportunity to merge two approaches in 
case of heavier nuclei. In the framework of the a-cluster 
model a nucleus is considered as a core (an a- cluster liq- 
uid drop, which grows with A) and a molecule [26l l27j . 
The last nucleon is supposed to be in the mean field of 
the molecule. So, with growing A the last nucleon po- 
tential stays unchanged, but the Coulomb potential for 
the last proton grows. In FIG. [3] a sketch of the spherical 
symmetrical potential is presented. In such a represen- 
tation the nuclear field of the nucleus comes from addi- 
tion of the potentials of the core and the molecule, so 
the requirement of larger internal part of the nucleus po- 
tential is fulfilled here. The parameters of the molecule 
potential can be found under the EPN condition for the 
nucleus corresponding to the molecule. The center of 
mass of the molecule is shifted towards the periphery on 
the value A, which grows with A. At < r < A the 
single particle wave function tf>(r) — and at r > A 
the wave function ip{r) — Cip'(r') where ip'(r') is the 
wave function of the last nucleon in cms of the molecule 
placed on the periphery of the Coulomb potential of a 
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FIG. 3: Local W-S potential with the Coulomb potential for 
the last proton in a nucleus with Z > 15 and he proton's wave 
function 



nucleus A and r' = r — A. The C is determined by the 
normalization of the functions J °° ip' \r') 2 r' 2 dr' = 1 and 
tp(r) 2 r 2 dr = 1. It will bring some other values of 
ANC and SF for heavier nuclei. Such representation will 
remove the mentioned above discrepancy and will provide 
a proper value Rq. Besides, the function will help in solv- 
ing the long standing problem of selecting optical poten- 
tial parameters used for description of input and output 
clastic channels of the reaction, because the internal part 
of the amplitude of the reaction will be naturally cut off 
at small radii, which will make differences of the values 
of the parameters of different optical potential sets less 
important for the calculated cross sections. Technically 
the potential can be used in numerical solutions of the 
Shrodinger equation without other difficulty. The param- 
eter A can be estimated in the framework of the a-cluster 
model or it can be found in the well-depth procedure at 
the other parameters fixed to adjust the experimental 
value of the single particle binding energy. The shell ef- 
fects in the binding energy, which are restricted within 
a few MeV, may come from varying the number of the 
a-clusters of the molecule from two to five. Practical 
application of this idea is a subject of another article. 

A simple phenomenological proof of this representation 
can be found in the values of the experimental binding 
energy of the last neutron in the symmetrical nuclei heav- 
ier than the nucleus 12 C. They group around two values 
e n =15 MeV for even Z and e n =11 MeV for odd Z 
( see Table II). So it means that the last neutron has 
two kinds of links in dependence on whether there is a 
single rtp-pair or not. The energy of neutron separation 



in odd(2')-odd (A) stable nuclei (according to the model 
these nuclei consist of a core made of a-clusters, nn-pairs 
placed in the core, and the molecule of 15 N) is also within 
few MeV around the value of the neutron separation en- 
ergy in the nucleus 15 N e n —11 MeV. For example, for 
the nuclei 19 F, 27 A1, 35 C1, 89 Y, 141 Pr and 209 Bi e„=10, 
13, 13, 12 and 9 MeV which also says about similar con- 
ditions of binding of the last neutron determined by the 
links with the single np-pair and three a-clusters in its 
close vicinity. 



V. CONCLUSION 

For some symmetrical nuclei A = B + n/p with 5< 
Z <14 the W-S potential parameters satisfying the EPN 
condition have been found. The Coulomb radius, which 
is the crucial parameter determining the other param- 
eters in such calculations, is obtained by using the last 
proton position radius R p in the cms of a nucleus A. The 
values R p are estimated from analysis of binding energies 
in the framework of the a-cluster model. The rms radii 
for the last proton < r 2 > 1 / 2 in case of the nuclei with 
Z <14 are in agreement with R p . In cases of pure pe- 
ripheral reactions with one nucleon transfer reactions like 
A(d, t)B or B( 3 He, d)A the error of spectroscopic factors 
obtained by means of the DWBA analysis with standard 
potential parameters is estimated. It is shown that the 
error for some nuclei is larger than 20%. 

Analysis shows that the W-S single particle bound 
state potential satisfying the EPN condition in case of 
the nuclei with Z >15 produces rms radius for the last 
proton bigger than R p . The radius of the potential in- 
creases with a leap and the potential radii are obviously 
too big to be consistent with the experimental radii. The 
standard parameters for these nuclei provide better fit- 
ting of rms radii to R p than the parameters of the EPN 
condition. But the potential radii are not consistent with 
the experimental charge radii. The better result is pro- 
vided by domination of the proton potential over neutron 
one, which does not seem reasonable. All this allows one 
to conclude that for the nuclei with Z > 15 the W-S po- 
tential is inappropriate to represent single particle bound 
state potential. 

According to the a-cluster model the potential well of 
the last nucleon is to be of a smaller radius due to the 
short range of nuclear force and it should be shifted to 
the nuclear periphery. An idea is proposed that for the 
DWBA analysis of one nucleon transfer reactions in case 
of the nuclei with Z > 15 the nuclear potential well for 
the last nucleon should be represented as a local W-S 
potential covering the periphery of the nucleus. 
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